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•  A  multi-scale  thermo  electrochemical  model  for  lithium  ion  cells  is  developed. 
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•  Predict  the  temperature  and  current  distribution  across  the  surface  of  the  cell. 

•  Compare  the  model  and  test  results  for  the  lithium  ion  cells. 
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Lithium  iron  phosphate  (LFP)  and  lithium  manganese  oxide  (LMO)  are  competitive  and  complementary 
to  each  other  as  cathode  materials  for  lithium-ion  batteries,  especially  for  use  in  electric  vehicles.  A  multi 
scale  multi-dimensional  physic-based  model  is  proposed  in  this  paper  to  study  the  thermal  behaviour  of 
the  two  lithium-ion  chemistries.  The  model  consists  of  two  sub  models,  a  one  dimensional  (ID)  elec¬ 
trochemical  sub  model  and  a  two  dimensional  (2D)  thermo-electric  sub  model,  which  are  coupled  and 
solved  concurrently.  The  ID  model  predicts  the  heat  generation  rate  (Q/i)  and  voltage  (I/)  of  the  battery 
cell  through  different  load  cycles.  The  2D  model  of  the  battery  cell  accounts  for  temperature  distribution 
and  current  distribution  across  the  surface  of  the  battery  cell.  The  two  cells  are  examined  experimentally 
through  90  h  load  cycles  including  high/low  charge/discharge  rates.  The  experimental  results  are 
compared  with  the  model  results  and  they  are  in  good  agreement.  The  presented  results  in  this  paper 
verify  the  cells  temperature  behaviour  at  different  operating  conditions  which  will  lead  to  the  design  of  a 
cost  effective  thermal  management  system  for  the  battery  pack. 

©  2014  Elsevier  B.V.  All  rights  reserved. 


1.  Introduction 

The  niche  market  of  the  battery  electric  vehicle  (BEV)  is  based 
on  reducing  the  production  cost.  In  order  to  reduce  the  cost  by  50% 
by  2030,  current  and  future  chemical  and  control  technologies 
must  be  successfully  incorporated.  Thermal  management  plays  a 
vital  role  in  the  pursuit  of  this  goal  1].  A  medium  range  BEV  with 
30  kWh  of  electric  energy  had  in  2011  a  battery  pack  with  a  mass  of 
300  kg.  This  is  expected  to  reduce  to  180  kg  by  2020  [2].  From  the 
battery  cost  perspective,  the  battery  cells  currently  cost 
$400  kWh-1  however  the  battery  pack  is  more  than  double  that 
cost  due  to  the  complex  thermal  and  electric  management  system 
required.  Research  and  Development  (R&D)  is  estimated  to  reduce 
the  cost  of  the  battery  pack  to  $300  kWh-1  by  2030  from  the 
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current  $1000  kWh-1  [3].  Most  of  the  current  R&D  is  on  improving 
the  energy  density  of  BEV  cells  via  new  chemistries,  rather  than 
investigating  thermal  aspect  of  the  cells.  The  development  of  new 
electrolytes,  cathode  and  anode  materials  is  the  core  of  cell  level 
studies  [4-6].  The  thermal  characteristic  behaviour  of  the  battery 
cell  is  also  investigated  in  different  papers  such  as  7].  The  majority 
of  them  studied  small  capacity  cells  (less  than  10  Ah  capacities)  and 
ignored  the  temperature  distribution  of  the  cell.  The  cell  di¬ 
mensions  in  BEV  application  are  more  than  10  times  those  of  the 
small  capacity  cells.  Therefore,  the  temperature  distribution  should 
be  considered  in  any  thermal  study. 

The  first  step  to  design  a  cost  effective  thermal  management 
system  for  the  battery  electric  vehicles  (BEV)  is  to  understand  the 
thermal  behaviour  of  the  battery  cell  as  the  main  source  of  heat. 
Having  reviewed  thermal  researches  on  the  lithium-ion  cell  applied 
in  BEV,  it  is  observed  that  researchers  studied  the  heat  generation 
and  heat  transfer  including  accumulation  and  dissipation  in 
different  levels.  Researchers  who  applied  the  first  principle 
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electrochemical  models  focused  on  the  heat  generation  terms  and 
left  the  heat  dissipation  part  for  the  battery  pack  level  study.  On  the 
other  hand,  researchers,  who  worked  on  the  heat  transfer  study  in 
batteries,  assumed  simply  a  resistive  heating  as  the  heat  source  in 
their  heat  transfer  model  and  focused  on  the  heat  dissipation  and 
temperature  distribution  through  battery  packs  [8].  However,  the 
thermal  issues  of  the  lithium-ion  cells  applied  in  EV  have  not  been 
fully  studied  concurrently.  This  paper  focuses  on  the  thermal 
behaviour  of  the  batteries  in  a  BEV  at  the  cell  level.  The  model  is  a 
first  principle  model  which  accounts  for  operating  temperature 
effect  on  the  performance  of  the  cell.  In  order  to  show  the  validity  of 
the  model  results,  the  lithium-ion  cells  are  experimentally  inves¬ 
tigated  and  the  results  are  compared  by  the  model  results  and  the 
model  results  are  in  good  agreement  with  the  experimental  mea¬ 
surements.  Both  aspects  of  the  thermal  study,  heat  generation  and 
dissipation,  are  considered  in  the  thermal  modelling.  The  heat 
generation  rate  is  calculated  and  the  model  is  capable  of  studying 
the  heat  dissipation  in  3D.  The  heat  generation  and  dissipation 
rates  are  dependent  variables  and  therefore  the  two  models  are 
coupled  to  detect  more  precisely  the  heat  generation  term  and 
temperature  distribution  throughout  the  cell.  Further,  the  cell 
component  properties  as  available  in  the  literature  are  function  of 
temperature  and  state  of  charge  and  are  incorporated  in  to  the 
model. 

2.  The  case  studies 

Lithium-ion  batteries  comprise  a  family  of  battery  chemistries 
that  employ  various  combinations  of  anode  and  cathode  materials. 
Lithium  iron  phosphate  (LFP)  and  lithium  manganese  oxide  (LMO) 
are  competitive  and  complementary  to  each  other  as  cathode 
materials  for  lithium-ion  batteries,  especially  for  use  in  electric 
vehicles.  LMO  batteries  have  been  used  in  electric  tools  and  ma¬ 
jority  of  on  road  electric  cars  such  as  Nissan  Leaf  use  this  type  of  the 
lithium-ion  chemistry;  LFP  battery  cell  also  uses  in  GM  electric 
vehicle,  Chevrolet  Volt.  However,  the  case  studies  in  this  paper  are 
not  representative  of  any  particular  electric  cars  and  their  typical 
data  are  given  in  Table  1.  The  case  study  battery  cells  have  the 
pouch  structure  and  are  dedicated  for  electric  vehicle  applications. 
The  fundamental  configuration  of  the  cells  is  shown  in  Fig.  1,  which 
is  not  to  scale  and  is  for  illustration  purposes.  The  principle  parts  of 
the  battery  cell  are  negative  electrode  collector  (Cu),  negative 
electrode,  separator,  positive  electrode  and  positive  electrode  col¬ 
lector  (Al).  This  configuration  repeats  several  times  (about  10 
times)  and  results  in  a  battery  thickness  of  approximately  7  mm  (Z- 
direction). 

3.  The  model  definition 


Fig.  1.  The  fundamental  configuration  of  the  battery  cell. 

been  applied  by  the  researchers  who  focus  on  the  cell  level  study. 
Therefore,  with  a  1 D  electrochemical  model,  most  of  the  reactions 
are  determined  without  undue  complexity  [9  .  This  approach  has 
proved  to  be  a  reliable  approach  for  projects  in  which  the  compo¬ 
nent  material  improvement  is  not  subject  of  the  study. 

The  proposed  model  consists  of  two  sub  models,  the  ID  elec¬ 
trochemical  sub  model  and  2D  thermo-electric  sub  model,  which 
are  coupled  and  solved  concurrently.  The  ID  model  can  predict  the 
heat  generation  rate  ( Qh )  and  voltage  (V)  of  the  battery  cell  through 
different  load  cycles.  The  2D  thermo  electric  model  of  the  battery 
cell  accounts  for  temperature  distribution  and  current  distribution 
across  the  surface  of  the  battery  cell.  The  model  from  a  numerical 
methods  perspective  is  a  finite  element  method  and  both  ID  and 
2D  models  are  discretised  into  computational  cells.  In  terms  of 
sequence,  the  model  input  is  a  load  cycle  which  is  representative  of 
the  expected  drawing  down  of  power/current  from  the  battery  cell 
with  time.  The  1 D  electrochemical  model  for  the  load  cycle,  solves 
the  mass  transfer,  charge  balance  and  electric  kinetic  equations 
concurrently  and  predicts  the  voltage  and  state  of  charge  for  each 
components  of  the  cell  (collectors,  electrodes  and  electrolyte).  The 
voltage  at  a  point  in  the  1 D  model  which  corresponds  to  the  tab  in 
the  battery  cell  is  then  mapped  to  the  tab  position  in  the  2D  thermo 
electric  model.  This  is  in  contrast  with  the  Gerver  [10]  combined 
model  which  is  composed  of  a  grid  of  ID  electrochemical  models 
coupled  to  a  thermal  model  and  a  2D  resistor  network  to  determine 
current  and  temperature  distributions  in  3D  dimensions.  However, 


The  electrochemical  reactions  in  the  battery  cell  are  defined  as 
heterogeneous  reactions  which  means  all  of  the  reactions  occur  at 
interfaces  between  electrode,  separator  and  electrode  in  Z  direc¬ 
tion,  see  Fig.  1.  On  the  other  hand,  the  high  electronic  conductivity 
of  the  current  collectors  (Al,  Cu)  bonded  to  the  electrodes  provide  a 
high  conductive  electric  layer  which  results  in  the  reaction  taking 
place  in  the  Z  direction.  Hence,  the  ID  electrochemical  model  has 


Table  1 

The  battery  cells  data. 


Parameters 

Lithium  iron 
phosphate  (LFP) 

Lithium  manganese 
oxide  (LMO) 

Capacity  (Ah) 

19.6 

20 

Voltage  (V) 

3.3 

4.2 

Cell  dimensions 

7  x  160  x  227 

8  x  177  x  240 

(mm) 

Fig.  2.  The  model  sequence. 
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based  on  the  input  from  ID  mode  that  are  the  tab  voltage  and  a 
volumetric  heat  generating  term,  the  2D  model  predicts  tempera¬ 
ture  and  current  for  each  computational  cells  and  then  weight 
average  of  the  temperature  across  the  surface  is  mapped  to  the  ID 
model  and  the  loop  is  repeated  till  the  solution  converges.  See 
Fig.  2. 

In  addition,  the  associated  total  heat  generation  rate  which  in¬ 
cludes  reversible  (reaction)  and  irreversible  (Ohmic  and  over  po¬ 
tential  terms),  maps  to  each  computational  cell  in  the  2D  thermo 
electric  model,  and  the  2D  model  solves  the  charge  balance  (Joule 
equation)  and  energy  equation  concurrently.  Hence  it  has  an  extra 
benefit  of  accounting  for  the  Ohmic  losses  across  the  surface  area  of 
the  cell  which  cannot  be  neglected  in  battery  cells  with  large  sur¬ 
face  area  such  as  EV  cells.  This  is  in  contrast  with  the  battery  model 
released  along  with  COMSOL  Multiphysics  4.2®  and  widely  was 
used  by  other  researchers  such  as  Xiao  and  Choe  [11  ]  in  which  only 
a  volumetric  heat  generation  term  is  applied. 

3.1.  The  ID  electrochemical  sub  model 


The  1 D  electrochemical  sub  model  is  based  on  Doyle  model  [12] 
which  has  been  built  by  COMSOLAa®.  The  model  was  developed  by 
White  et  al.  to  include  thermal  effects  [1  ].  Bernadi  et  al.  applied  the 
model  to  analyse  the  relaxation  behaviour  in  cells  designed  for 
hybrid  vehicle  propulsion  [13].  It  solves  six  coupled,  nonlinear 
equations  to  calculate  the  cell  voltage  and  current  as  a  function  of 
time,  and  to  calculate  the  potentials  in  the  electrolyte  and  electrode 
phases,  salt  concentration,  intercalated  lithium  concentration,  re¬ 
action  rate,  and  current  density  in  the  electrolyte  as  functions  of 
time  and  position  within  the  cell. 

The  charge  balance  0\  in  the  solid  phase  is  governed  by  the 
Ohm’s  law. 


v(-/<ffV0,)  =  ,~SaJloc  (1) 

k\  is  the  electronic  conductivity,  j\oc  is  the  local  charge  transfer 
current  density  (positive  for  anodic  reactions)  (A  m-2),  Sa  is  the 
specific  surface  area  (m2  m“3).  The  ionic  charge  balance  for  the 
liquid  phase  expresses  as: 


Vi 


aln/ 


9ln  c2 


[1  -  t+]V(ln  c2)  }  =  S;JU 


(2) 


The  subscripts  1  and  2  denote  to  the  solid  and  liquid  phases  of 
active  materials.  In  the  Doyle  model,  the  mass  balance  of  lithium 
ions  in  an  intercalation  particle  of  electrode  active  material  is 
described  by  Fields  law.  This  is  expressed  as  follows  by  introducing 
the  dimensionless  radius,  y  =  r/rp,  where  rp  is  the  particle  radius,  Di 
solid  phase  diffusivity,  and  c\  is  the  concentration  of  lithium  in  the 
active  material. 


(3) 


The  material  balance  for  the  binary  electrolyte  in  the  liquid 
phase  is  given  by: 


dc? 


+  V<!  “  -  T  0  -  4)  \  =  0 


(4) 


C2  is  the  cation  concentration  in  electrolyte,  e  denotes  the  porosity 
(volume  fraction  of  electrolyte),  and  Dlff  the  effective  diffusivity  in 
the  electrolyte.  Butler-Volmer  equation  describes  the  electrode 
kinetics. 


Jloc  —  *0 


(5) 


where 


V  =  <t>l  -^2-£ref(Cl,  surf) 


(6) 


i o  =  k< 


:0 \jc2  (c1. 


cl,  surf 


h 


surf 


(7) 


z'o  denotes  the  exchange  current  density,  and  Eref  refers  to  the 
electrode  equilibrium  potential  of  particle,  which  is  a  function  of 
the  concentration  at  the  particle  surface,  c\urf.  In  addition,  /c0  de¬ 
notes  the  reaction  rate  constant,  and  Ciimax  is  the  maximum  sur¬ 
face  concentration.  The  activation  over  potential,  r /,  for  all 
electrode  reactions  in  the  electrode  receives  an  extra  potential 
contribution. 


3.2.  Temperature  dependency  of  parameters 

One  of  the  accepted  formulae  for  temperature  dependence  of 
reaction  rate  is  the  Arrhenius  equation.  [1,14]  The  Arrhenius 
equation  is  best  seen  as  an  empirical  relationship  and  is  used  to 
model  the  temperature  variation  of  diffusion  coefficients  and 
many  other  parameters  [15].  However,  measured  correlations  for 
some  parameters  also  proposed  [16].  The  lithium-ion 
cell  parameters,  temperature  and  concentration  dependence 
equations  of  electrolyte  conductivity,  salt  diffusivity  and  ther¬ 
modynamic  factor  are  taken  from  Ref.  [17].  The  diffusivity  of 
negative  electrode  is  temperature  dependent  and  is  given  by  Ref. 
[7]. 

•  Electrolyte  conductivity  (k\): 


k,  =  O.OOOld  (  -  10.5  +  0.668  x  10“3d  +  0.494*10-6c? 
+0.074T-  1.78  x  10“5d  x  T-8.86  x  10-10c?  x  T 

-  6.96  xlO“5T2  + 2.8  x  10“8Ci  xl2)2  (8) 

•  Electrolyte  salt  diffusivity: 

A  -4.43-f - ^ ^  )-0.22xl0-3cZ  , 

D,  =  10~4  X  10  4-229-5x,0-3d,/  (9) 

•  Activity  dependence: 


din/ 
©In  cl 


(0.601  -  0.24(l0-3d)°5  +  0.982(1 
-  0.0052(7  -  294))  (lO“3d)  ’  5) 


(10) 


•  Negative  electrode  intercalation  diffusivity: 


DSn  =  1.4523  x  10 


13  exPr  68025  7  (  1  1 

F\8. 314472  \318  T 


(ID 


Table  2 

The  2D  model  parameters. 
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Z  (pm) 

X  (mm) 

Y  (mm) 

Density 

(g  cm-3) 

Heat  capacity 

Usk-1) 

Thermal  conductivity 
(W  mK1) 

Electric  conductivity 
(S  ITT1) 

Positive  collector 

20 

177 

240 

2.7 

0.9 

238 

3.83E+07 

Positive  electrode 

150 

177 

240 

1.5 

0.7 

5 

139 

Negative  collector 

10 

177 

240 

8.96 

0.385 

398 

6.33E+07 

Negative  electrode 

145 

177 

240 

2.5 

0.7 

5 

100 

Separator 

52 

177 

240 

1.2 

0.7 

1 

J[D OD,T)  <  1.5 

Equivalent  parameters 

Thermal  conductivity: 

2.14 

0.678 

27.23 

J[  DOD.T) 

In  X,Y  direction  parallel 

Ay 

27.23 

In  Z  direction  series 

Az 

3.39 

3.3.  Energy  equation 

The  thermodynamics  of  the  cell  is  based  on  the  thermal  energy 
balance  over  an  elementary  volume  in  a  battery.  The  differential 
equation  that  describes  temperature  distribution  in  the  battery 
takes  following  conservation  form 

QT  n 

pCp0F  =  Avr  +  (?g-<?d  (12) 


through  each  sheet  in  series,  while  heat  flowing  along  the  stack 
flows  through  each  layer  in  parallel.  Because  the  thickness  of  each 
layer  is  very  small,  significant  computational  time  is  consumed 
when  attempting  to  resolve  heat  flowing  in  each  layer.  Instead,  the 
battery  unit  cell  is  typically  modelled  as  a  uniform  layer  with 
anisotropic  thermal  conductivity  and  uniform  heat  generation  as 
follows: 

Thermal  conductivity  in  X  and  Y  direction  (/;  denotes  the 
thickness  of  layer  i): 


qd  is  the  rate  of  heat  exchanges  between  the  battery  cell  and 
surroundings. 


,  EM 

XT 


(16) 


<Jd 


^cell 

^cell 


(J  —  Ts) 


(13) 


h  is  the  heat  transfer  coefficient  (W  nrr2  K-1),  T  is  the  cell  tem¬ 
perature  and  Ts  is  the  surrounding  temperature.  Aceii  and  VCeii  are 
the  surface  area  and  volume  of  the  battery  cell. 

The  heat  source  term  (qg)  in  the  battery  cell  is  defined  by  Ref. 

[18] 


%  —  ^2  asjlnj  (0s  —  0e  —  Uj)  +  ^2  asjinjT  -^jr  +  CeffV0sV0s 
j  j 

+  /<effV0eV0e  +  kffv In  C\  V0e 

(14) 

The  first  term  on  the  right  hand  side  is  potential  deviation  in 
the  control  volume  from  the  equilibrium  potential  (irreversible 
heat).  The  second  term  represents  the  entropic  effects  (reversible 
heat).  The  third  term  shows  the  Ohmic  heat  in  the  solid  phase  and 
the  last  two  terms  represents  the  Ohmic  heats  in  the  solution 
phase. 


3.4.  The  2D  thermo-electro  model 


The  2D  model  is  an  equivalent  model  representing  the  3D  ge¬ 
ometry  of  the  battery  cell  and  the  thermal  and  electrical  properties 
of  the  material  are  used  to  calculate  the  equivalent  properties.  The 
2D  mode  parameters  are  given  in  Table  2.  The  battery  properties  in 
the  X  and  Y  directions  are  considered  to  be  in  parallel  and  in  the  Z 
direction  is  series.  See  Fig.  1  for  direction.  The  battery  is  a  layered 
cell  containing  multiple  materials  with  different  thicknesses  and 
thermal  properties.  See  Fig.  1.  As  a  result,  the  thermal  conductivity 
of  the  battery  is  anisotropic.  The  composite  charge  conductivity  is 
calculated  as  follows 


k  = 


RT  RT  \  1 
io,anF  +  i0cnFJ 


(15) 


The  heat  flows  perpendicular  to  the  X  direction  and  flows 


Thermal  conductivity  in  z  direction 

**  =  E  i  07) 

The  heat  capacity  for  the  equivalent  cell  is  volume  averaged  as 
follows 


pCp  =  5>C P)ih/J2h 


(18) 


3.5.  Thermal  boundary  condition 

The  battery  cells  are  studied  to  determine  the  cells  behaviour  in 
the  absence  of  any  active  cooling  system.  The  boundary  condition 
has  a  significant  effect  on  the  temperature  distribution  of  the  cell. 
The  experiment  and  the  model  in  this  project  are  set  up  in  a  same 
way.  The  ambient  temperature  is  measured  in  the  test  and  used  as 
the  surrounding  temperature  (Ts)  in  the  model,  equation  (13).  The 
cell  during  the  test  was  placed  vertically  in  a  condition  without 
having  any  conductive  heat  transfer  with  the  surrounding.  In  the 
model,  a  heat  transfer  coefficient,  h,  is  considered.  The  heat  transfer 
coefficients  for  all  surfaces  of  the  cell  are  considered  as  a  constant 
number  equal  to  6  W  m-2  K-1  [1,19].  This  number  is  obtained  by 
calculating  the  Nusselt  number  of  the  buoyancy  flow  around  the 
cell.  However,  the  model  is  capable  to  set  up  different  heat  transfer 
boundary  conditions  and  is  used  to  study  different  heat  transfer 
coefficients  and  ambient  conditions. 

4.  The  experimental  measurements 

The  high  capacity  battery  cells,  20  Ah,  at  different  load  cycles  are 
examined.  Three  battery  cells  with  the  same  specification  are  tested 
concurrently  and  shown  that  the  test  results  are  repeatable.  To 
ensure  the  validity  of  the  test  results,  the  uncertainty  and  level  of 
confidence  of  the  test  results  are  calculated  and  compared  to  the 
equipment  accuracy.  Three  fundamental  measurements,  tempera¬ 
ture,  terminal  voltage  and  current  throughout  the  load  cycle  are 
measured.  The  parameters  are  typically  measured  at  periodic 
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intervals  during  a  battery  test  sequence.  The  15  thermocouples  are 
set  as  a  grid  across  the  surface  of  the  cells  to  measure  the  temper¬ 
ature.  The  tests  are  performed  at  ambient  temperature  which  is  also 
monitored  throughout  the  tests.  In  order  to  have  repeatable  results, 
three  battery  cells  with  the  same  specification,  chemistry  and  con¬ 
struction  are  tested  concurrently.  In  this  project,  data  acquisition 
intervals  of  4  s  were  applied  and  the  sample  data  reviewed  carefully 
to  be  sure  that  all  samples  are  independent.’  See  Fig.  3. 

Inverse  Distance  Weighting  (IDW)  method  [20]  is  used  to  plot 
temperature  contours  from  thermocouples  data.  The  IDW  method 
estimates  the  values  of  an  attribute  by  using  a  linear  combination  of 
values  at  sampled  points  weighted  by  an  inverse  function  of  the 


measurement  data  especially  in  high  current  charge/discharge 
rates  (bigger  than  1C).  For  example,  the  cell  after  a  10  min  discharge 
with  5C  (100  A)  reaches  a  temperature  of  45  °C  and  the  corre¬ 
sponding  temperature  for  the  model  in  this  condition  is  45.5  °C.  In 
addition,  the  temperature  drops  and  reaches  the  ambient  temper¬ 
ature  at  the  same  time  in  both  the  model  and  the  test  data. 

There  is  a  subtle  difference  between  the  predicted  temperature 
and  the  experimental  data  during  the  charge/discharge  process 
with  the  rate  (less  than  1C),  see  Fig.  5.  For  example,  the  model 
predicts  temperature  of  28  °C  for  1C  discharge  rate  which  is  4° 
lower  than  that  measured  in  the  test.  Flowever,  the  difference  is 
minimum  at  the  peak  point  and  the  model  reaches  the  ambient 


Fig.  3.  The  test  set  up. 


distance  from  the  point  of  interest  to  the  sampled  points.  In  addi¬ 
tion,  a  thermal  camera  is  applied  to  detect  hot  spots  and  the  tem¬ 
perature  distribution  on  the  surface  of  the  battery  cells.  The  two  set 
of  experiment  results,  thermocouple  data  and  the  thermal  image, 
show  the  same  gradient  across  the  sample’s  surface  and  identified 
the  hot  zone  on  the  sample. 

5.  Results  and  discussion 

The  cells  are  tested  with  a  load  cycle  lasting  90  h  with  range  of 
constant  current  charge  and  discharge  rates  from  0.2C  (0.2  times 
nominal  current  that  is  20  A)  to  10C.  See  Fig.  4.  The  cell  voltage  and 
current  were  acquired  for  each  time  step,  4  s.  Temperature  over  the 
surface  of  the  cell  was  measured  by  15  Thermocouple  probes  with 
their  distribution  on  the  cell  surface  shown  in  Fig.  3. 

The  same  load  cycle  is  applied  for  the  model  and  test.  The  model 
prediction  for  average  temperature  of  the  LFP  cell  surface  and  the 
measured  average  temperature  are  illustrated  in  Fig.  5.  In  the  model 
it  is  assumed  that  convection  was  from  the  large  plane  cell  surfaces 
with  a  heat  transfer  coefficient  of  h  =  6  W  m-2  K  1  [19],  and  an 
ambient  temperature  which  is  constant  and  equal  to  20  °C.  Fig.  5 
shows  the  model  prediction  is  in  good  agreement  with  the 


Fig.  4.  The  LFP  test/model  load  cycle. 


temperature  at  the  same  time  and  this  is  important  as  it  shows  that 
the  model  thermal  mass  compensates  for  the  difference  in  pre¬ 
diction  of  the  peak  value.  The  LFP  cell  study  shows  that  the  pre¬ 
dicted  temperature  in  low  charge/discharge  rate  (less  than  1C)  are 
not  as  good  as  high  rate  in  agreement  with  the  experimental  data. 
The  difference  between  heat  generation  terms  in  low  and  high 
current  rate  could  contribute  to  make  the  difference  between  the 
predicted  and  measured  temperature.  In  high  current  rate,  the 
irreversible  term  dominates  the  other  heat  generation  terms  and 
the  Ohmic  term  in  the  model  correctly  represents  the  Ohmic  heat  in 
the  battery  cell.  In  the  low  current,  the  reversible  term  which  is  the 
sum  of  producible  reversible  work  and  entropic  heating  from  the 
reaction  contributes  more  on  the  generated  heat.  The  model’s 
predicted  reversible  heat  is  not  the  right  representative  of  revers¬ 
ible  heat  which  occurred  during  the  measurement.  This  hypothesis 
is  verified  by  conducting  another  set  of  experimental  measurement 
with  different  lithium-ion  chemistry.  Therefore,  a  lithium  manga¬ 
nese  oxide  battery  (LMO)  cell  is  selected  for  further  study  which  is 
explained  in  the  flowing  section. 

The  LMO  battery  cell  which  is  chosen  for  the  study  has  the  same 
capacity  as  the  LFP  cell,  20  Ah.  Fig.  6  shows  the  model  prediction  for 
average  temperature  of  the  LMO  cell  surface  and  the  measurement 
of  average  temperature.  In  the  model  it  is  assumed  that  convection 
was  from  the  large  plane  cell  surfaces  with  a  heat  transfer  coeffi¬ 
cient  of  h  =  6  W  m~2  K-1  [19],  and  an  ambient  temperature  which  is 
considered  by  the  measured  ambient  temperature  throughout  the 
test.  As  it  can  be  seen  in  Fig.  7,  the  model  voltage  prediction  is  in 
good  agreement  with  the  measurement  either  high  or  low  current 
discharge  rates.  The  electrochemical  model  can  provide  the  polar¬ 
isation  curve  in  the  accuracy  of  95%.  In  addition,  temperature 
dependent  parameters  and  the  heat  generation  and  dissipation 
terms  considered  in  the  model  predict  the  average  temperature 
with  the  accuracy  93%.  The  measured  temperature  and  predicted 
rise  considerably  over  the  discharge  period  with  the  same  trend. 
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Fig.  5.  Comparison  between  mean  temperature  of  the  model  result  (LFP  cell)  and  the  test. 


Fig.  6.  Comparison  between  mean  temperature  of  the  model  result  (LMO  cell)  and  the  test. 


Fig.  7.  Voltage  curves  for  the  LMO  cell  at  different  discharge  rates  (1C  =  20  Ah). 


Fig.  8.  Mean  temperature  variation  with  depth  of  discharge  for  the  LMO  cell  at 
different  discharge  rates  (1C  =  20  Ah). 


However,  the  measured  temperature  jump  and  fall  5%  when  the  cell 
20%  is  discharged  which  is  not  seen  in  the  model  especially  for  high 
currents  rates.  See  Fig.  8. 

The  results  demonstrate  that  at  low  current  rates  (0.2C),  the 
LMO  cell  performance  is  better  in  terms  of  maintaining  the  voltage 
during  discharge.  These  low  current  rates  would  correspond  to  cell 
use  in  EV  whilst  the  high  discharge  rates  are  associated  with  use  in 
hybrid  vehicles,  and  hence  a  short  life. 

In  order  to  obtain  the  temperature  distribution  across  the  cell 
surface,  thermocouples  data  and  thermal  camera  are  used.  The 
predicted  temperature  contours  across  the  surface  of  the  cell  are 
compared  with  the  thermal  camera  image  at  3C  discharge;  Fig.  9.  It 
shows  that  both,  thermal  camera  and  the  model  prediction,  deter¬ 
mine  the  same  pattern  for  the  temperature  distribution  across  the 
surface  of  the  cell.  Overall  the  temperature  difference  across  the  cell 
at  the  end  of  discharge  is  about  2  °C  and  this  corresponds  to  drawing 
60  A  current  from  the  cell  for  20  min.  The  maximum  temperature 
occurs  at  the  top  of  the  cell  near  the  tabs  and  the  contours  are 
symmetric.  However,  the  temperature  contour  depends  on  how  set 
up  the  thermal  boundary  condition.  In  the  presented  results,  it  is 
assumed  that  the  heat  dissipates  from  the  edges  and  surfaces  which 
is  an  appropriate  estimation  for  natural  convection  based  on 
buoyancy  flow  and  calculating  Nusslet  number  [19].  Temperature 
and  voltage  distributions  are  analogues  [21  ],  and  thus  the  contours 
on  Fig.  8  represent  voltage  distribution  across  the  cell  as  well. 

6.  Conclusion 

A  multi  scales  physics  based  model  is  developed  to  study  the 
battery  cell  behaviour  and  the  model  results  are  compared  by  the 
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Fig.  9.  Comparison  between  temperature  contour  of  the  model  result  and  the  test. 


experimental  results.  The  model  voltage  prediction  is  in  agreement 
with  the  measurement  either  high  or  low  current  discharge  rates. 
The  electrochemical  model  can  provide  the  polarisation  curve  in 
the  accuracy  of  95%.  In  addition,  the  model  predicts  the  average 
temperature  with  the  accuracy  93%. 

Two  different  lithium  chemistries  are  applied  and  the  results 
show  the  model  can  predict  the  temperature  and  voltage  of  lithium 
polymer  cell  (LMO)  more  accurate  than  for  the  lithium  iron  phos¬ 
phate  (LFP)  cell.  The  study  results  show  the  LFP  cell  parameters’ 
dependency  functions  should  be  obtained  experimentally  or  the 
functions  should  be  adapted  to  apply  for  high  capacity  cells  as  most 
of  them  have  not  been  applied  for  EV  application. 

For  BEV  which  operates  at  a  nominal  low  current  rate,  the  order  of 
magnitude  of  heat  source  through  the  cell  is  less  than  850  W  m-3 
during  discharge  process ;  this  amount  of  heat  can  be  managed  to  keep 
the  cell  at  desired  operating  temperature.  BEV  battery  is  more  likely 
encounter  problems  to  keep  it  warm  at  the  low  end  of  the  tempera¬ 
ture  range,  than  to  require  cooling  at  high  ambient  temperature.  The 
temperature  difference  across  the  cell  at  the  end  of  discharge  is  about 
2  °C  and  this  corresponds  to  drawing  60  A  current  from  the  cell  for 
20  min  in  ambient  temperature  of  20  °C.  The  maximum  temperature 
occurs  at  the  top  of  the  cell  near  the  tabs  and  the  contour  pattern 
depends  on  the  boundaries’  heat  transfer  condition. 

It  is  recommended  to  study  the  battery  cell  thermal  behaviour  in 
typical  and  realistic  driving  cycles,  continue  the  concurrent  study  of 
heat  generation  and  dissipation  in  EV  batteries  pack  level  and 
compare  performance  of  thermal  management  techniques  to 
maintain  the  cell  in  the  desired  operating  temperature  in  terms  of 
lightweight  and  less  sophisticated  cooling  systems. 


(7eff  effective  matrix  conductivity  (1  (Qm)"1) 

0i  potential  in  the  matrix  phase  (V) 

02  potential  in  the  solution  phase  (V) 
k  electronic  conductivity  (S  m-1) 

Sa  specific  surface  area  (m2  m-3) 

R  gas  constant  (8.3143  J  mol-1  I<-1) 

y  fraction  of  lithium  in  the  positive  electrode 

7  Bruggeman  coefficient 

lc  thicknesses  of  cathode  electrode  (m) 

T  temperature  (I<) 

F  Faraday’s  constant  (96487  C  mol-1) 

/  ionic  activity  coefficient 

t+  transport  number 

e  porosity  (volume  fraction  of  electrolyte) 

Dfff  effective  diffusivity  in  the  electrolyte  (m2  s-1) 
io  exchange  current  density  (A  m-2) 

Eref  electrode  particle’s  equilibrium  potential  (V) 
Cit  SUrf  particle  surface  concentration  (mol  m-3) 
ko  reaction  rate  constant 

Ci, max  maximum  surface  concentration  (mol  m-3) 

7]  activation  overpotential  (V) 

Di  electrolyte  salt  diffusivity  (m2  s-1) 

din  f /din  cl  activity  dependence 
Dsn  Negative  electrode  intercalation 

diffusivity  (m2  s-1) 

Uref  open  potential  of  the  electrode  (V) 

x  fraction  of  lithium  in  the 

negative  electrode 

la  Thicknesses  of  anode  electrode  (m) 


List  of  symbols 

Ci  concentration  of  the  electrolyte  (mol  m-3) 

cs  surface  concentration  in  spherical  particle  (mol  m-3) 

Cp  specific  heat  (J  kg-1  K-1) 

jioc  transfer  current  density  (A  m-2) 

q  volumetric  heat  generation  rate  (J  m-3  s-1) 

L/j  equilibrium  potential  (V) 

keff  effective  solution  conductivity  (1  (Qm)_1) 

X  thermal  conductivity  (W  m-1  K-1) 

p  density  (kg  m-3) 


Appendix  1.  The  electrochemical  model  equations 

kf  =  ki^  ,  i  =  1,2  (A-l) 

Df  =  Di  £(  ,i=l,2  (A-2) 

From  Ref.  [22] 

Urefn  =  -0. 1 6  + 1 .32  exp(  -  3x)  + 1 0  exp(  -  2000x)  (A-3 ) 

From  Ref.  [23] 
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^ref,LMO 


-23.53y5  +  52.43y4  -  44.86y3  +  18.4y2  -  3.35y 
+  3.16 

(A-4) 


From  Ref.  [16] 


=  344.1347148  exp(-32.9633287x  +  8.316711484)/(l 

+  749.0756003  exp(-34.7909964x  + 8.887143624)) 
-0.8520278805x  +  0.36229929x2 +  0.2698001697 

(A-5) 


^  =  4.131274309  exp(0.571536523y) 

di/LMO  V 

+  1.281681 122  sin(- 4.991 6739y) 

-  0.090453431  sin(  -  20.9669665y  + 12.5788250) 

-  .0313472974  sin(31 .7663338y  -  22.4295664) 

-  4.14532933  +  8.1471 13434y  -  26.064581y2 
+  12.7660158y3  -  0.184274863 

x  exp(  -  (y-  0.51 694351 68)/0.04628266783)2^ 

(A-6) 

From  Ref.  [10] 


^ref,LFP 


-5.76  x  10+  +  2.07  x  104y3  -  2.79*104  xy2 
+  1.67  x  104y  -  3.77  x  103 

(A-7) 


=  — 0.35376y8  +  1 ,3902y7  -  2.2585y6  +  1 ,9635y5 

/  LFP 

-  0.9871 6y4  +  0.28857y3  -  0.046272y2 
+  0.0032 158y  -  1.9186  x  1CT5 

(A-8) 
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